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^ ■ This paper investigates Magnetic Reynolds number dependence of the "self-similar 

O ' evolution model" (Nitta et al. 2001) of fast magnetic reconnection. I focused my attention 

| 7j \ on the flow structure inside and around the reconnection outflow, which is essential to 

determine the entire reconnection system (Nitta et al. 2002). The outflow is consist of 
several regions divided by discontinuities, e.g., shocks, and it can be treated by a shock- 
tube approximation (Nitta 2004). By solving the junction conditions (e.g., Rankine- 
Hugoniot condition), the structure of the reconnection outflow is obtained. Magnetic 
reconnection in most astrophysical problems is characterized by a huge dynamic range of 
its expansion (~ 10^ for typical solar flares) in a free space which is free from any influence 
of external circumstances. Such evolution results in a spontaneous self-similar expansion 
lO which is controlled by two intrinsic parameters: the plasma-/? and the magnetic Reynolds 

O I number. The plasma-/? dependence had been investigated in our previous paper. This 

paper newly clarifies the relation between the reconnection rate and the inflow structure 
^! just outside the Petschek-like slow shock: As the magnetic Reynolds number increases, 

P \ strongly converging inflow toward the Petschek-like slow shock forms, and it significantly 

reduces the reconnection rate. 

a • Subject headings: Earth — MHD — Sun: flares — ISM: magnetic fields 



^1 1. INTRODUCTION of the coupling process between different scales, 

H : Magnetic reconnection is one of probable candi- j-^- ^^^^ ^he microscopic scale process re- 

■ ■ ■ ' dates for very powerful energy conversion mech- latmg to the origm of anomalous resistivity (~ 

anism in magnetized plasma systems, e.g., solar 1° ["^] ^olar corona) govern the macroscopic 

flares and geomagnetospheric substorms. Major- ^^fl^ P^of^^^ m which enormous magnetic energy 

ity believes that magnetic reconnection is a uni- ^61^^^^^ (~ 1° H ^yPi^^l ^^^6 l^^P) ^^^^ 

versal mechanism also for recent energetic topics, important, but still an open question. One may 

e.g., accretion disk X-ray emission of young stellar ^^'""^ ^^^^ ^^"o^s ^^a* extremely small diffu- 

objects (YSOs; see Koyama et al. 1994; Hayashi ^^o^ ^^g/o^ go^^™^ the huge entire reconnection 

et al. 1999) and galactic ridge X-ray emissions system (e.g., idea of the ' fractal diffusion region ' ; 

(GRXEs; see Koyama et al. 1986; Tanuma et al. ^^^^ ^hibata & Tanuma 2001).^ The aim of the se- 
1999). 



ries of our papers is to clarify the relationship and 



There are so many phenomenological studies causality between these two different scales m a 

of magnetic reconnection mainly by numerical regime of MHD (see subsection 4.4). In this work, 

simulation researchers, however, essential open our attention is focused on the magnetic Reynolds 

questions still remain, not only regarding the number dependence of the flow structure and the 

microscopic physics of the anomalous resistivity, reconnection rate. 

but also the macroscopic magnetohydrodynami- ^e consider a magnetic energy conversion m 

cal (MHD) structure. Especially, understanding antiparallel magnetic equilibrium configura- 

tion called a "current sheet system. ' In such 
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the case, the electric resistivity plays an impor- 
tant role: It breaks the frozen-in and leads to the 
magnetic diffusion. The magnetic diffusion can 
release the magnetic energy by the Ohmic dissi- 
pation. However since the Spitzer resistivity of 
the most cases of astrophysical plasma systems is 
considerably small, the magnetic diffusion speed 
should be very small. This means that the mag- 
netic diffusion itself is not efficient for a quick en- 
ergy conversion like in the solar flares. We here 
introduce an essential non-dimensional parameter 
so-called the "magnetic Reynolds number Rem 
which is defined as the ratio of the Alfven wave 
transit timescale across the system to the mag- 
netic diffusion timescale. The most of astrophysi- 
cal plasma is highly conductive (e.g.. Re m ^ 10 
for typical solar corona), thus the energy conver- 
sion speed by the Ohmic dissipation is very slow 
(the energy conversion timescale is of the order 
of 10^[yr]!!: This is too long to explain the solar 
flares) . 

It is important that, even in such the case, the 
processes involving magnetic reconnection can re- 
lease the magnetic energy at a very large power. 
The mostly successful reconnection model for as- 
trophysical applications is the Petschek model 
(see Petschek 1964). Although the reconnection 
of magnetic field lines itself is caused by the non- 
uniform magnetic diffusion, the main energy con- 
version is owing to the propagation of the X- 
shaped slow shocks in the Petschek model. Con- 
sequently, the Petschek model can convert the 
magnetic energy in a very short time (the wave 
propagation timescale estimated as of the order 
of 10^ [sec] for typical solar flares) almost indepen- 
dent of the magnetic Reynolds number. We call 
such a very powerful process in which the mag- 
netic diffusion timescale does not determine the 
entire energy conversion as the fast reconnection. 
As discussed above, the Petschek model might 
be the most significant model to explain actual 
reconnection phenomena, but the author thinks 
there is an essential lack to be a real astrophysi- 
cal reconnection model. The following part of the 
introduction is focused on this point. 

We must note that many cases of the actual 
magnetic reconnection in astrophysical systems 
usually grow over a huge dynamic range in its spa- 
tial dimension. For example, the initial scale of 
the reconnection system can be defined by the ini- 
tial current sheet thickness, but this is too small 
to be observed in typical solar flares. We do not 



have any convincing estimate of the scale, but if 
we estimate it to be of the order of the ion Lar- 
mor radius, it is extremely small (~ 10° [m] in 
the solar corona). Finally, the reconnection sys- 
tem develops to a scale of the order of the initial 
curvature radius of the magnetic field lines (~ 10^ 
[m] ~ 1.5% of the solar radius for typical solar 
flares). The dynamic range of the spatial scale is 
obviously huge (~ 10'' for solar flares). For geo- 
magnetospheric substorms, their dynamic range 
of growth is also large (~ 10^ for substorms). 
Such a very wide dynamic range of growth sug- 
gests that the evolution of the magnetic reconnec- 
tion should be treated as a development in "free 
space", and that the external circumstance does 
not affect the evolutionary process of a magnetic 
reconnection, at least at the expanding stage just 
after the onset of reconnection. 

Even a system is not completely free from the 
influence of the external circumstances, we can 
approximately treat it as a spontaneously evolv- 
ing system if the evolution timescale (which is es- 
timated as the Alfven transit time for the flnal 
system scale because the spontaneous expansion 
speed of the reconnection system is equivalent to 
the fast-mode propagation speed, see Nitta et al. 
2001; hereafter paper 1) is much smaller than the 
timescale imposed by the external circumstances 
(e.g., the convection timescale). For typical sub- 
storms, the evolution timescale of the reconnec- 
tion is of the order of second or minute while 
the convection timescale to compress the current 
sheet is of the order of hour. In such cases, the 
external circumstances simply play a role of trig- 
gering the onset of a reconnection and the evo- 
lution itself is approximately free from the influ- 
ence of the external circumstances. Though, of 
course, the exceptional cases in which the exter- 
nal circumstances intrinsically influence the evo- 
lution can arise for the particular situation (e.g., 
the cases that very fast convection drives the re- 
connection), the author believes that it is worth 
to establish a new model applicable to the cases 
which are free from the external influences. 

However, no reconnection model evolving in a 
free space has been studied. We should note that 
the most of previous theoretical (e.g., Petschek 
1964, Vasyhunas 1975, Priest & Forbes 1986) and 
numerical works on a reconnection treated it as a 
boundary problem strongly influenced by external 
circumstances. Our interest is focused on a spon- 
taneous evolutionary process of a magnetic recon- 
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nection in a free space as a macroscopic instabil- 
ity in a current sheet system. This evolutionary 
reconncction model can explain a coupling pro- 
cess between different scales, i.e., a microscopic 
scale represented by the ion Larmor radius and 
a macroscopic scale represented by the scale of 
external circumstances. 

The author and collaborators had been studied 
such the problem of two-dimensional MHD recon- 
nection in a series of papers. A non-linear evo- 
lution of the spontaneous reconncction in a free 
space is numerically simulated (two-dimensional 
MHD simulations with artificially enhanced lo- 
calized resistivity) in paper 1. The result of our 
numerical simulation clearly shows a self-similar 
expansion of the reconnection system. Owing to 
technical reasons (restriction of CPU time, mem- 
ory etc.), we could confirm the stable expansion of 
the reconnection system only in a dynamic range 
of 10^ — 10^ in spacial scale of the expansion. This 
dynamic range is obviously insufficient because 
actual dynamic range seems to be 10^ — 10^ as 
discussed above. In order to supplement the lack 
of the dynamic range, we studied the inflow region 
of the reconnection system semianalytically by 
so-called the "Grad-Shafranov approach" (Nitta 
et al. 2002; hereafter paper 2) with a boundary 
condition along the Petschek-like X-shaped slow 
shock (Nitta 2004; hereafter paper 3). We ob- 
tained a self-similarly expanding solution which 
is fairly consistent with our numerical result in 
paper 1. Therefore we can conclude that the self- 
similar solution of an expanding reconnection sys- 
tem is stable with any ideal MHD modes and will 
continue to expand until the system scale becomes 
to be comparable to the typical scale of the exter- 
nal environment, e.g., the diameter of a magnetic 
flux loop. In paper 3, we studied a structure of the 
reconnection outflow which provides the bound- 
ary condition along the Petschek-hke slow shock 
by a shock-tube approximation. This boundary 
condition is very important because it determines 
the structure of the inflow region (see paper 2). 

The spontaneously expanding reconnection sys- 
tem has only two parameters: the plasma-/? value 
and the magnetic Reynolds number. We could 
clarify the plasma-/? dependence for a very wide 
dynamic range of plasma-/? in paper 3 (a brief 
summary is in the next section). This enables us 
to apply our self-similar reconnection model to a 
wide variation of phenomena which may have dif- 
ferent plasma-/?. In this paper, we study the mag- 



netic Reynolds number dependence of the recon- 
nection system. While majority believes that a 
spatially localized anomalous resistivity is impor- 
tant to realize a fast reconnection (e.g., Yokoyama 
& Shibata 1994), we do not have any established 
model for the mechanism and a convincing esti- 
mated value of the anomalous resistivity (Coppi 
& Friedland 1971, Ji ct al. 1998. Shinohara et 
al. 2001) yet. Thus, it is important to study 
the magnetic Reynolds number dependence of the 
properties of the reconnection system. We have 
obtained the following results. 

We define the effective magnetic Reynolds num- 
ber as i?*^ = Vao/V^*/ where Vao is the Alfven 
speed at the asymptotic region (the region far 
from the current sheet) and V^*j is the actual 
magnetic diffusion speed at the diffusion region. 

1) Strongly converging inflow 

As i?*^ increases, the convergence of the inflow 
toward the Petschek-like slow shock@becomes sig- 
niflcant. This leads that the inflow velocity and 
the magnetic field line at the infiow region tend 
to be parallel to each other. 

We here define the reconnection rate R* as the 
non-dimensional reconnection electric filed (the 
ratio to the product of Vao and Bq [the magnetic 
field strength at the asymptotic uniform region]). 

2) Reconnection rate 

As a result of the converging infiow, the recon- 
nection electric filed decreases as i?*^ increases. 
Finally we find R* (x 1/ Rl^ at the large magnetic 
Reynolds number limit. 

This paper is organized as follows. We state 
a model of the reconnection outfiow as a kind of 
shock-tube problem in section 2. The basic equa- 
tions for the outfiow structure and the numerical 
procedure to solve the equations are listed also 
in this section. By solving the basic equations, 
we obtain the fiow structure also the reconnec- 
tion rate. The results are shown in section 3. We 
discuss the properties of these results in section 
4. 

2. NUMERICAL SCHEME 

2.1. Concept 

As we mentioned in paper 3, the reconnection 
outfiow in our spontaneous model has a struc- 
ture involving several discontinuities like a shock- 
tube problem. We solve the reconnection out- 
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flow region by approximating the problem as a 
kind of shock-tube problem (see details in the 
next subsection). In paper 3, we assumed a non- 
converging inflow toward the slow shock approx- 
imating the original Petschek model, and dis- 
cussed the structure of the outflow region as a 
function of the ambient plasma-/? value. Conse- 
quently, we obtain the plasma-/? dependence of 
the spontaneous inhalation speed (see figure 1) 
and the reconnection rate (see figure 2). These 
figures are essentially the same with the corespon- 
dent figures in paper 3, but newly plotted by a 
more precise calculation. 

In this paper we are willing to study the mag- 
netic Reynolds number dependence of the outflow 
structure by the above shock-tube approximation 
in the regime of ideal (non- resistive) MHD. One 
may think that it is very curious because we can- 
not treat the magnetic Reynolds number in the 
regime of ideal MHD as in this work. We should 
note, however, that if we can evaluate the recon- 
nection electric field Ez (note that Ez is uniform 
in the vicinity of the reconnection point from the 
Ampere's law because this region is almost sta- 
tionary in our self-similar model) and the mag- 
netic field B* around the diffusion region from 
the shock-tube problem, we can estimate the in- 
fiow speed v* (= E^/B*) toward the diffusion re- 
gion. This inflow speed is in balance with the 
magnetic diffusion speed v^if at the diffusion re- 
gion. Providing that the Alfven speed Vao at 
the asymptotic region to be flxed as a normal- 
ization value as in paper 3, inverse of the inflow 
speed toward the diffusion region shows the ef- 
fective magnetic Reynolds number i?*^^ (-Rem = 
^Ao/v^if = Vao/v*). Thus we can treat the mag- 
netic Reynolds number in a context of ideal MHD. 
We would stress again that our attention is fo- 
cused on the magnetic Reynolds number depen- 
dence of the reconnection system, i.e., the outflow 
structure and the reconnection rate. 

2.2. Model 

We present a schematic picture of the reconnec- 
tion outflow in our self-similar evolution model 
(see figure 3). The coordinates are defined as fol- 
lows: a;-axis is parallel to both the current sheet 
and initial anti-parallel magnetic field, y-axis is 
perpendicular to the current sheet, and 2;-axis is 
parallel to the current sheet, but is perpendicu- 
lar to the initial magnetic field (hence, dz=0 in 
this two-dimensional problem). The essence of 



the model is quite the same as the model in pa- 
per 3. A similar structure of the reconnection 
outflow is shown in figure 9 of Ugai (1999) as a 
simulation result. Because of symmetry with re- 
spect to X— and |/— axes, we treat only the region 
X, y > 0. The outflow is composed of two differ- 
ent plasmas which have different origins. These 
two plasmas arc touch in contact at the point 
X = Xc (the contact discontinuity). A plasmoid 
forms around the contact discontinuity. The rear- 
half region {x < x^. the reconnection jet) is filled 
with the reconnected plasma coming from outside 
the current sheet (the inflow region). The front- 
half region {x > Xc- the plasmoid) is flUed with 
the original current sheet plasma. 

The entire outflow is surrounded by a slow shock 
which has a complicated 'crab-hand' shape (see 
Abe & Hoshino 2001). The Petschek-hke figure-X 
shaped slow shock (an oblique shock) is elongated 
from the diffusion region with a slight opening an- 
gle 9. There is a reverse fast shock (an almost 
perpendicular shock) inside the reconnection jet 
{x = Xf). In front of the plasmoid, a forward V- 
shaped slow shock (an oblique shock) forms. The 
opening angle and the locus (the crossing point 
with X— axis) are (j) and Xg-, respectively. 

The entire structure including the several dis- 
continuities is analogous to the one dimensional 
shock-tube problem. The forward shock (a V- 
shaped slow shock) and the reverse shock (a fast 
shock) are formed by the collision of the recon- 
nection jet and the original current sheet plasma, 
and propagate in both directions. Between these 
two shocks, the contact discontinuity forms. We 
approximate this reconnection outflow as a quasi 
one-dimensional problem in order to solve it an- 
alytically. Such an approximation may be valid 
near the axis, because the system is symmet- 
ric with respect to the a;— axis. 

We focus our attention on the quasi-one- 
dimensional problem along the inflow stream line 
and the reconnection outflow. We treat the figure 
L-shaped region x,y = finite 3> if the inflow 
stream line is parallel to the y— axis where D is 
the initial current sheet thickness. Note that this 
region apparently tends to coincide with the x— 
and I/— axes in the self-similar stage which is a 
very late stage from the onset when we observe 
the evolution in a zoom-out coordinate [see pa- 
per 1]. If the inflow stream lines are inclined 
owing to the convergence of the inflow, we con- 
sider the modifled L-shaped region along the in- 



S. Nitta 



5 



flow stream line (see SSL of figure 7) in stead of 
the y—sods. Each region between two neighboring 
discontinuities is approximated to be uniform. We 
note about the up-stream region p just above the 
X-shaped slow shock that the region between the 
slow shock and the separatrix field line (X-shaped 
field line reaching the X-point) is also approxi- 
mated to be uniform. In this region, each recon- 
nected field line has an almost straight shape and 
crosses the X-slow shock while each field line has 
hyperbolic shape in the region above the separa- 
trix field line. In paper 3, we presumably suppose 
a non-converging inflow in the region p for the 
spontaneous inhalation, but if inhalation speed 
is reduced when i?*^ increases, it may result in 
an alteration of the inflow structure. Hence, we 
should take a possibility of convergence or diver- 
gence of the inflow in the region p into our prob- 
lem (sec below and subsection 2.4). 

The region 3 between the contact discontinuity 
and the V-shaped slow shock looks to be non uni- 
form. Since we do not know how we can solve the 
two dimensional structure of this region analyti- 
cally, we approximate this region to be uniform. 

The quantities denoting the initial uniform 
equilibrium at the asymptotic region are gas 
pressure Pq, mass density po and magnetic field 
strength Bq. The plasma-/? value at the asymp- 
totic region is defined as Pq = Po/{Bq/[2ijl]) where 
fi is the magnetic permeability of vacuum. In the 
rest of this paper, we use the normalization of 
physical quantities as like in paper 2. We de- 
fine units for each dimension as follows: (unit of 
velocity)= Vao = Bo/^/JIp^ (Alfven speed at the 
asymptotic region), (unit of the length) = VAot 
where t is the time from the onset of reconnec- 
tion, (unit of mass density) = po, (unit of magnetic 
field)= Bq, (unit of pressure)= Po/2 ■ PoV^o- 

We set the quasi-one-dimensional shock-tube- 
like problem as follows. The system has 22 un- 
known quantities: Pp, pp, v^p, B^p, Byp, 6, Pi, pi, 

Vi, Bi, Xf, P2, P2, V2, B2, P3, Ps, Vy3, ^xS, ByS, 

(f) and Xs where P*, p*, v*, B^ denote the pres- 
sure, density, velocity, magnetic field, respectively 
(Note Xc = Vx3 = V2 because no mass flux passes 
through the contact discontinuity). The suffices 
p, 1, 2 or 3 denote the region divided by the dis- 
continuities. The suffices x or y denote the vector 
components. Quantities 9, Xf, (f) and Xg denote 
the inclinations of the X-shaped slow shock, the 
locus of the fast shock, the inclination of the V- 
shaped slow shock and the locus of the V-shaped 



slow shock (crossing point with a;— axis), respec- 
tively. 

These unknowns should be related to each other 
via conditions coming from the integrated form 
of conservation laws (i.e., the Rankine-Hugoniot 
[R-H] conditions) or other relations. The set of 
relations is listed in the next subsection. 

There is an essential difference from the model 
in paper 3. In paper 3, we assumed (or imposed) 
the transversal component of the inflow velocity 
vanishes [v^p = 0) and solved the vertical com- 
ponent Vyp of the inflow velocity (and other un- 
knowns) . We here treat artiflcially con- 
trollable variable and its value is given by hand. 
The change of Vyp may induce a change of the 
strength of the fast-mode rarefaction. Since the 
inflow is induced by the fast-mode rarefaction, 
this change may result in a change of v^p- Hence, 
we must solve Vxp as a new unknown quantity. 

2.3. Basic equations 

There is no essential difference of the basic equa- 
tions for the outflow structure between the equa- 
tions used in paper 3 except that we must pre- 
pare the equations denoting the oblique inflow 
case Vxp 7^ 0. In the most cases, we must simply 
add several terms involving v^p (e.g., in equations 
3, 5, 6, 7, 13 and 22 in paper 3). 

In equation 1 (the frozen-in condition at the re- 
gion p), a physical insight to modify it is required. 
In the oblique {v^p 7^ 0) inflow case, we should 
note that the frozen-in condition requires that the 
ratio of the perpendicular component B± of the 
magnetic fleld with respect to the inflow velocity 
versus the mass density holds to be a conservative 
quantity. If we can presumably obtain the stream 
line conflguration all over the inflow region, we 
can state exact modified equation for the frozen- 
in condition, however, we never know it in the 
context of this quasi-one-dimensional shock-tube 
approximation. Hence we need an approximating 
model of the infiow stream line configuration: we 
here simply approximate it to be straight. Thus 
we obtain equation Al in appendix. 

According to the above model of the reconnec- 
tion outfiow, we obtain the following 22 equations 
for 22 unknown quantities (Detailed forms of each 
equation are listed in the appendix). 

l)relations between pre-X-slow shock and asymp- 
totic region 

Prozen-in condition [eql] 
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Polytropic relation [eq2] 

2) R-H jump conditions at X-slow shock 
Pressure jump [eq3] 

Density jump [eq4] 

Velocity jump (parallel comp.) [eq5] 

Velocity jump (perpendicular comp.) [eq6] 

Magnetic field jump (parallel comp.) [eq7] 

Magnetic field jump (perpendicular comp.) [eq8] 

3) R-H jump conditions at reverse fast shock 
Pressure jump [eq9] 

Density jump [eqf 0] 
Velocity jump [eqf 1] 
Magnetic field jump [eqf 2] 

4) Magnetic flux conservation at X-point [eqf 3] 

5) Force balance at contact discontinuity [eqf 4] 

6) R-H jump conditions at forward V-slow shock 
Pressure jump [eqf 5] 

Density jump [eqf 6] 
Velocity jump (parallel comp.) [eqf 7] 
Velocity jump (perpendicular comp.) [eqf 8] 
Magnetic field jump (parallel comp.) [eqf 9] 
Magnetic field jump (perpendicular comp.) [eq20] 

7) Boundary Condition at the tip of the outflow 
[eq2f] 

8) Magnetic flux conservation all over the outflow 
[eq22] 

We can solve the majority of these equations by 
hand, and by substituting the solutions into other 
equations, we can reduce equations. Finally, ten 
equations f, 3, 4, 5, 6, 7, ff, f4, 2f and 22 re- 
main as complicated nonlinear coupled equations 
for ten unknowns Pp, v^p, B^p, Byp, 9, Pi, pi, Xf, 
(f) and Xs- 

2.4. Numerical procedure 

We solve these coupled ten equations by an it- 
erative method (the Newton- Raphson method). 
In order to find the well converging solution of 
the nonlinear coupled equations, a precise initial 
guess of the unknowns is required. In general, 
this step to find an appropriate initial guess is 
the core of difficulty (this should be comparable 
to a "treasure hunting in ten-dimensional space" 
with an incomplete treasure map). Fortunately, 



this most difficult step had already been cleared 
in paper 3 for the case v^p — 0. 

We can start from the solution of the case 
Vxp = for an arbitrary value of (3o. Our most in- 
terest is focused on low-/3 cases typically appear in 
astrophysical applications. We demonstrate the 
numerical procedure for the case = O.Of as 
the reference case bccaiisc this is a typical value 
in solar corona and gcomagnetosphere. In the 
work of paper 3, we have already had the solu- 
tion in the spontaneous inhalation case {v^p — 0) 
for po = O.Of. 

First, we start with the solution of the case 
Vyp = -Vinsp{0.01) = -0.06398 (where Vinsp{0.01) 
is the spontaneous inhalation speed for Pq — 
O.Of) as the initial guess for the case ^yp — 
— ^jnsp(O.Of) + Avy, where Avy is the increment 
of Vyp. Once we find a converged solution of the 
Newton-Raphson procedure for the case Vyp = 
—Vinsp{0.01) + Avy, we treat it as the initial guess 
for the case Vyp — —Vinsp{0.01) + 2Avy, and we 
have successively obtained the solutions for differ- 
ent Vyp. Of course, Avy should be a small enough 
value to keep good convergence of the Newton- 
Raphson method. From several trials, we care- 
fully adopt the following two cases of the incre- 
ment: Avy = 10^^ and fO"'''. We have checked 
that the results of these two cases fairly coincide 
with each other. 

We obtain the result that, as \vyp\ decreases, the 
strength of the reverse fast shock decreases. At 
a critical value of \vyp\, the reverse fast shock no 
longer forms (the density jump ratio and the pres- 
sure jump ratio reduce to unity). At this point 
we exchange the numerical code to another one 
for the case of no reverse fast shock forms (see 
the last part of appendix), and we can continue 
the calculation. Above procedure is valid for any 
value of (3o. 

3. RESULT 

We have investigated the response of the recon- 
nection outflow structure to a variation of Vyp (y- 
component of the inflow velocity at the pre-slow 
shock region [region p] normalized by the ambient 
Alfven speed Vao). Usually, it shows the reconnec- 
tion rate if the transverse component B^p of the 
magnetic field is almost same as its ambient value 
Bq and x-component v^p of the inflow velocity is 
negligible {\vxp\ <^ bj/p|)- We must note, however, 
that if a signiflcant converging inflow arises or the 
transverse magnetic component reduces from its 



S. Nitta 



7 



ambient value, we need an alternative estimation 
of the reconnection rate. Actuary a significant 
converging inflow and decrease of the transverse 
magnetic component realize when \vyp\ decreases 
(see figure 4). Thus we introduce the effective 
reconnection rate 

which denotes the reconnection electric field 
normalized by VaqBq. 

The magnetic Reynolds number i?em is defined 
as the ratio of the ambient Alfven speed to the 
magnetic diffusion speed. We should note that the 
shock-tube approximation adopted here is within 
a regime of ideal MHD, so we cannot treat the 
diffusion region itself in this scheme. Although 
one may think that we cannot argue the dif- 
fusion speed in this scheme at all, actually we 
can estimate it by a simple approximation as fol- 
lows. In the vicinity of the reconnection point, 
the system is approximately in a steady state in 
the self-similarly expanding solution. This fact 
leads, from the Ampere's law, the electric field 
E* = —VypBxp + VxpByp is approximately uniform 
around the diffusion region. The diffusion speed 
w^jj must be balanced with the local infiow speed 
v*^ in the stationary state, and is estimated 
as v*^ = E*/B* = {-VypBxp + VxpByp)/B*, where 
B* is the magnetic field just outside the diffusion 
region. We suppose B* ~ B^p in order to approx- 
imately satisfy the magnetic flux conservation. 

Thus, we obtain the effective diffusion speed 

V*dif = V-n ~ {-VypBxp + VxpByp)/ Bxp , (2) 

and define the effective magnetic Reynolds 
number 



We select the |/-component Vyp of the infiow 
velocity toward the slow shock as the artificially 
controllable variable for our numerical parameter 
survey. Starting at the solution with its sponta- 
neous inhalation speed Vyp = —VinspiPo) toward 
the slow shock and gradually decrease its absolute 
value with a finite decrement (sec section 2.4), and 
numerically solve the coupled equations listed in 
section 2.3 by the Newton-Raphson method for 
each value of Vyp. Thus we obtain a series of solu- 
tions for the outfiow structure. Note that i?*^ is 
a monotonically increasing function of Vyp in our 
result (see figure 4), so we adopt i?*^ as the con- 
trol parameter in the following discussion. The 
result is as follows; 

1) Converging infiow 

The converging speed \vxp\ (note Vxp < in 
our result which denotes a converging inflow) 
of the inflow drastically increases as the mag- 
netic Reynolds number R*^ increases (see figure 
4). The directions of the infiow velocity and the 
magnetic field tend to be parallel as the converg- 
ing speed increases (see figure 6). 

2) Magnetic Reynolds number dependence of re- 
connection rate 

As i?*^ increases, first the reconnection rate R* 
increases owing to the increase of Bxp (see figure 
4), and then R* decreases owing to that the direc- 
tions of the infiow velocity and the magnetic field 
tend to be parallel (see figure 4). This causes a 
reduction of the reconnection electric field. Thus, 
the reconnection rate R* decreases as i?* oc 1/ -R*„j 
because Bxp ~ 1 in the large R^^^ region (see fig- 
ure 4). 



Km = "^hlif = Bxp/ {-VypBxp + VxpByp) . (3) 

Note Vao = 1 in our normalization. Hence, we 
obtain 

R* = Bxp/ Rl^ ■ (4) 

It is very important to clarify the relation be- 
tween R* and R^^. For this purpose, we treat 
i?*^ as a control parameter instead of Vyp. The 
effective magnetic Reynolds number (i^em) depen- 
dence of the spontaneously expanding reconnec- 
tion system in a free-space is investigated for the 
cases Po — 0.01 and Po — 0.2. 



4. DISCUSSION 

4.1. Channel flow structure in the inflow 

For large magnetic Reynolds number (say. 
Rim — 20), a significantly converging infiow to- 
ward the slow shock arises. It can grow to be 
almost parallel to the magnetic field lines as R*„^ 
increases. In such the case, the infiow toward the 
diffusion region must be very slow because the 
reconnection rate (note R* = Vp x Bp) is very 
small, even though the Pctschck-likc slow shock 
still strongly intakes the plasma in order to sat- 
isfy the Rankine-Hugoniot condition. This means 
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that the inflow forms a channel flow structure: rel- 
atively large speed converging inflow toward the 
slow shock and very slow inflow toward the diffu- 
sion region (see figure 7). 

It is possible to observationally estimate the 
reconnection rate from the inflow Alfven Mach 
number (see Yokoyama et al. 2001, Isobc ct al. 
2005). If we can obtain the Alfven speed and the 
inflow speed perpendicular to the inflow magnetic 
fleld by observation (i.e., if we observe a motion 
of a bright fllament-like structure toward the X- 
point, supposing it shows a fleld line motion, we 
can estimate the normal speed of the fleld fine mo- 
tion) , we can estimate the reconnection rate as a 
good approximation. We must note, however, the 
inflow Alfven Mach number itself does not denote 
the reconnection rate if the inflow is inclined with 
respect to the magnetic fleld lines. Such the sit- 
uation may occur when the flow is signiflcantly 
converging (see section 3). In such the case, the 
reconnection rate might be rather small even the 
inflow speed is very large. We must note that, in 
order to estimate the reconnection rate, we need 
not only the inflow Alfven Mach number, but also 
the direction of the magnetic fleld relative to the 
inflow velocity. 

4.2. The role of the diffusion region 

In a very short period (of the order of the Alfven 
transit time over the current sheet thickness) just 
after the onset of the reconnection, the evolution 
strongly depends on the resistivity model. Such 
information propagates by MHD waves. In the 
self-similar stage, since the effect of the resistivity 
model remains only in a thin layer in the vicin- 
ity of the FRWF, it does not influence the entire 
structure of our self-similar solution. When the 
induced inflow sufficiently evolves and Petschek- 
like slow shock forms, the slow shock begins to 
play major roles for spontaneous structure forma- 
tion. The region around the diffusion region with 
a finite spatial scale seems also to be determined 
by the slow shock and tends to be stationary state 
as time proceeds. The resistive region plays only 
passive roles in this stage. These dynamics had 
been partially clarified by our numerical simula- 
tion (sec Nitta et al. 2001). The discussion in 
this paper aims to clarify the spontaneously de- 
termined structure lead by the shock-tube prob- 
lem in the region described by ideal MHD. The 
following is just a conjecture in the present sta- 
tus, but I believe that it can be clarified in our 



near future works. 

In the self-similar stage, the size of the diffu- 
sion region (similar to the current sheet thickness 
or the ion Larmor radius) is infinitesimally small 
comparing with the entire system scale. Hence, 
the most of the inflow does not pass through the 
diffusion region, but pass through the Petschek- 
likc slow shock to form the reconnection outflow. 
In such situation, it may be trivial that the entire 
structure should be determined not by the diffu- 
sion region, but by the Petschek-like slow shock. 
The only role of the diffusion region (but is not 
less meaningful) is to restrict the diffusion speed 
of the magnetic fleld lines. This effect of the diffu- 
sion region is clarifled by this paper (see sections 
3 and 4.1). 

In my conjecture, the structure of the diffu- 
sion region in the self-similar stage is passively 
determined in a circumstance of the surround- 
ing non-resistive region which is determined by 
the shock-tube problem involving the Petschek- 
like slow shock as is treated in this paper. In this 
case, the environmental circumstance which is de- 
termined by the shock-tube problem will play a 
role of a kind of boundary condition to determine 
the diffusion region. If we focus our attention into 
the diffusion region, this situation is very similar 
to the case of the "driven reconnection" which is 
also determined by the external boundary condi- 
tion. 

This intuition should be certifled by solving the 
diffusion region under some matching condition 
with the environmental ideal MHD region. We 
must realize, however, that the structure of the 
diffusion region cannot be described by MHD but 
by some particle scheme, because the scale of the 
intended diffusion region is of the order of the 
ion Larmor radius. This is obviously out of the 
limit of MHD approximation. To inquire further 
into the matter would lead us into that special- 
ized area of the particle reconnection, and such 
a digression would undoubtedly obscure the out- 
line of our argument in the regime of the MHD 
reconnection. 

4.3. Petschek-type slow shock formation 

In the shock-tube approximation adopted in 
this paper, we assume the steady formation of 
the Petschek-type figure-X shaped slow shock as 
a presupposition of the problem since we are in- 
terested in the fast reconnection regime. 
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One may think this is curious because the en- 
tire system is strongly time-dependent. Though 
wc cannot find any definite answer in the present 
status, the author would show the following con- 
jecture. First, we must realize that the central re- 
gion with a finite spatial scale around the diffusion 
region is in a steady state in the self-similar stage 
even the entire system is self-similarly expanding. 
In this case, the central region will be controlled 
by the surrounding stationary circumstance which 
is determined by the shock-tube problem (see the 
previous subsection). In naive intuition, we may 
expect a resultant stationary diffusion region and 
steady slow shock formation in the central region. 
The further study of this problem lies outside the 
scope of this paper as discussed in the previous 
subsection. 

We must note that the assumption of steady 
slow shock formation might not be appropriate 
for the case with a very small reconnection rate. 
The converging speed \vxp\ drastically increases 
as the magnetic Reynolds number increases 
(see figure 4), and it exceeds unity at i?*^ ~ 50. 
This is very curious. As long as the infiow is 
spontaneously induced by the reconnection sys- 
tem, its speed may not exceed unity for the case 
of low /^o because the Alfven speed is the maxi- 
mum speed when magnetic energy is completely 
converted to the bulk kinetic energy. This curi- 
ous result might owing to our presupposition that 
the Petschek type slow shock always forms: In or- 
der to keep to form the slow shock, the Rankine- 
Hugoniot condition requires unrealistic very high 
speed infiow. 

In actual reconnection, if the resistivity is very 
small and hence the magnetic Reynolds number 
becomes large (say, R*^ > 50), the Petschek type 
slow shock may not form. Thus, reconnection 
may change its type to another one with no slow 
shock (e.g., the Sweet-Parker type; Sweet [1958], 
Parker [1963]) or X-O-X type one (see subsection 
4.6). For the Sweet-Parker type of reconnection, 
we cannot adopt a shock-tube approximation, and 
it should be beyond the scope of this paper. 

4.4. Coupling between microscopic and 
macroscopic processes 

The size of the diffusion region seems to be ex- 
tremely small if wc estimate it by of the order 
of ion Larmor radius (~ 10*^ [m] for typical solar 
corona). From a naive intuition, it is very curious 
that such a very small diffusion region governs the 



huge entire reconnection system (e.g., the radius 
of magnetic flux tube, say ~ 10'^[m] for typical 
solar flares). Shibata & Tanuma (2001) proposed 
a new idea of "fractal reconnection" that effective 
diffusion region has a fractal structure in which 
total scale of the diffusion region is not so small, 
and is composed of spatially self-similar small 
structures. Owe self-similar evolution model, how- 
ever, gives an alternative picture about this prob- 
lem as follows. 

The self-similar reconnection model grows in a 
very wide dynamic range of its spatial scale. It 
may be triggered when the system scale is of the 
order of ion Larmor radius which is the scale of 
the current sheet thickness to drive an anoma- 
lous resistivity. The self-similar solution can ap- 
proximate the actual phenomenon when the re- 
connection system scale is much larger than ion 
Larmor radius (say, > 10^[m]). This self-similar 
solution can describe the system until the expand- 
ing system scale reaches an environmental proper 
scale which is imposed by external circumstances. 
Hence, we must note that the self-similar recon- 
nection model is describing the coupling process 
between the microscopic scale and the macro- 
scopic scale. 

Needless to say, the major energy release is de- 
termined by the largest structure of the recon- 
nection system. There are two different points 
of view about how is the largest structure de- 
termined. In the most of previous works, it 
is treated to be governed by external circum- 
stances through the boundary condition (the ex- 
ternally "driven" reconnection, see Vasyliunas 
1975, Priest & Forbes 1986, Sato & Hayashi 1979, 
etc.) . In this scheme, the external circumstances 
directly determine the large scale structure inde- 
pendent of the microscopic process of the diffu- 
sion region. In another scheme, the large scale 
structure is determined by intrinsic causes which 
properly included in the system itself, so we call 
it the "spontaneous" reconnection. Ugai & Tsuda 
(1977) first showed a spontaneous model: the re- 
sistivity model denoting microscopic processes of 
the anomalous resistivity determines the entire 
system. Our self-similar evolution model is an ad- 
vanced spontaneous reconnection model extended 
to the time-dependent reconnection which may be 
actual in astrophysical applications. 

The evolution just after the onset of reconnec- 
tion will depend on detail properties of the re- 
sistivity model, i.e., size and shape of the resis- 
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tive region, value and distribution of the resistiv- 
ity inside the diffusion region, etc. Anyway, once 
reconnection starts, its size expands in the fast- 
mode rarefaction wave (FRW) propagation speed. 
When the system scale (i.e., the scale of the wave 
front FRWF of FRW) becomes much larger than 
the initial scale (the initial current sheet thick- 
ness ~ ion Larmor radius), the system tends to 
be described by our self-similar solution. After 
this stage, the reconnection will proceed unlim- 
itedly in a self-similar way. The self-similar solu- 
tion depends only on two parameters (the plasma- 
(3 value at the asymptotic region and the ef- 
fective magnetic Reynolds number i?*^), and is 
independent of other properties. For example, 
effects of above mentioned detailed properties of 
the resistivity model is contracted within a very 
narrow region in the vicinity of FRWF. As the 
self-similar expansion proceeds, the influence of 
detailed properties is quickly diminished. Be- 
cause the self-similar solution well approximates 
the system as long as its scale is much larger than 
the initial current sheet thickness (~ 10*^ [m]), our 
self-similar solution is valid after the system scale 
becomes, say 10^ [m]. This is very small compar- 
ing with the final system scale (~ 10'^[m]). Af- 
ter that stage, the structure does not change no 
more but simply expands self-similarly. Thus, we 
should realize that the self-similar solution which 
is determined when the system scale is very small 
(~ 10^[m]) governs much larger structure in a 
later stage by the manner of self-similar expan- 
sion. Even the system scale reaches some exter- 
nal scale, the properties of the self-similar solution 
hold for a very long period (~ several hundred 
times the fast-mode transit time, see section 4.5 
of paper 1). 

4.5. Steady vs. time- dependent and driven vs. 
spontaneous 

In "time-dependent" reconnection models start- 
ing from an equilibrium state, a plasmoid forms 
around the contact discontinuity between the 
reconnection jet and the initial current sheet 
plasma. This plasmoid is a result of the inter- 
action between these two different plasmas. In 
our "spontaneous" model in a free-space, the en- 
tire system is self-consistently determined from 
this interaction (see Nitta et al. 2002 and Nitta 
2004). Thus, we do not need any external bound- 
ary condition in order to determine the resultant 
reconnection system. 



On the other hand, in "steady" models (see Va- 
syhunas 1975 or Priest and Forbes 1986), the in- 
formation of this interaction is completely lost. 
Thus, we need to impose the external boundary 
condition in order to determine the reconnection 
system. This is so-called a "driven" type recon- 
nection. A similar discussion about the difference 
between steady and time- dependent reconnection 
is in Ugai and Zheng (2005). 

Needless to say, the answer to the question 
which type is relevant depends on the circum- 
stance of the problem. The author had been 
discussed the relation between the spontaneous 
model and the driven model in paper 3 (see sec- 
tion 5.6) in detail. I here emphasize again that 
these two different models are not opposite sub- 
jects, but models for different phases of the entire 
evolutionary process of the astrophysical recon- 
nection. In my opinion, the "driven phase" fol- 
lows the "spontaneous phase" in usual astrophys- 
ical problems as follows. 

A large scale plasma convection (e.g., convec- 
tion of fiux tubes in the solar plasma) will form 
a current sheet system, and stores an enormous 
magnetic energy into the current sheet system. 
In usual cases of convection with the speed at 
very small Alfven Mach number, it will need a 
very long time comparing with the Alfven transit 
timescale for the energy storage (e.g., in typical 
solar fiarcs, it needs several days ~ 10^ [s]). As 
the result of current sheet thinning, once a well- 
localized anomalous resistivity is switched on in 
the current sheet, it triggers the onset of a fast 
reconnection. The evolution quickly moves to the 
self-similar expansion. The reconnection system 
expands at the speed of the FRW propagation 
which is estimated as the ambient Alfven speed 
(much faster than the convection speed). Hence, 
the evolution just after the onset is "spontaneous" 
one, because the timescale of the expansion is 
much smaller (~ 10^ [s] for typical solar fiares) 
than the timescale imposed by the external cir- 
cumstance (~ 10^ [s]). When the FRWF reaches 
the ambient circumstance, the evolution will start 
to be infiuenced by the external boundary con- 
dition, but the central energy conversion region 
holds as the spontaneous state during almost one 
hundred times the Alfven transit timescale (see 
section 4.5 of Nitta et al. 2001). After the central 
region is altered by the external boundary condi- 
tion, the system moves to the "driven" phase at 
last. 



S. Nitta 



11 



4.6. Possibility of another new type of solution 

Many people believe that the Petschek type fast 
reconnection can realize under the existence of the 
locally enhanced anomalous resistivity (Biskamp 
1986, Scholer 1989, Yokoyama & Shitaba 1994). 
We have not been able to identify the mechanism 
of the anomalous resistivity yet, but some kind of 
current driven microscopic instability should be 
the origin of the anomalous resistivity. In naive 
picture, such current driven anomalous resistivity 
may arise at the locus where the current sheet is 
strongly compressed. 

In this paper, we considered only the case that 
the diffusion region is fixed around the origin of 
the coordinate. If we apart from this presumable 
assumption, we may realize another possibility. 
From the discussion of the channel flow structure 
(see subsection 4.1), we can speculate the follow- 
ing new type of solutions. In the channel flow, the 
inflow speed toward the slow shock is much larger 
than it toward the diffusion region. This non- 
uniform inflow leads a non-uniform compression 
of the field reversal region. Since the compression 
outside the diffusion region is much stronger ow- 
ing to a strong inhalation by the slow shock, the 
reconnection point may move outside the present 
diffusion region. A bipolar reconnection outflow 
will be ejected from the new X-point. One out- 
flow towards the direction apart from the previous 
diffusion region, another towards the previous dif- 
fusion region. Thus, an 0-point will form around 
the original diffusion region. Consequently, a typ- 
ical X-O-X structure will form. We must note 
that this X-O-X structure is not a result of the 
tearing-mode instability but of the non-uniform 
inflow speed lead by the Petschek-like slow shock. 

The existence of this new solution shows a pos- 
sibility of the bifurcation of the solution. If the 
magnetic Reynolds number Rem is small enough 
(say, < 20), the Petschek-like X-point structure 
will be stable because no channel flow structure 
forms. When Rem increases and exceeds a critical 
value (~ several x 10), a remarkable channel flow 



structure will take place, and the solution will bi- 
furcate into two different types of solutions: the 
single X-point solution and the X-O-X solution. 

A similar discussion about non- uniqueness of 
the reconnection solutions was shown in Hameiri 
(1979). He apphed his model to flare phenomena 
as abrupt jumps to different branches of the so- 
lution. The answer to the question, in our case, 
which solution is stable will depend on the re- 
sistivity model (artificially fixed type around the 
origin or a current driven type) and the "history" 
of the time variation of the magnetic Reynolds 
number. It may result in a hysteresis behavior. 
An indication of the transition of the solution to 
the X-O-X structure is appeared in our prelim- 
inary numerical simulation. The details will be 
discussed in our following papers. 

We must realize, however, that we cannot com- 
pletely avoid these kinds of ambiguity arising from 
the non-uniqueness of the resistivity model as 
long as our study is restricted within the regime 
of MHDs because the resistivity cannot be de- 
termined from the MHD equations. Unfortu- 
nately, we have not had any convincing micro- 
scopic model of anomalous resistivity. I believe 
that these facts do not reduce the importance to 
study the MHD reconnection; on the contrary, we 
must study MHD reconnection more as an ele- 
mentary process in astrophysical plasma before 
applying to complicated global problems because 
we have not established complete astrophysical 
MHD reconnection model yet. 
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APPENDIX 

EQUATIONS FOR STRUCTURE OF THE RECONNECTION OUTFLOW 
The structure of the reconnection outflow is determined by the following equations. 



We assume that the region between the asymptotic region and the pre-shock region is filled with non- 
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resistive plasma. Hence, the magnetic flux is frozen into the induced inflow. We also assume a polytropic 
variation in the induced inflow because there is no violent process in the fast-mode rarefaction. Thus, 
we impose 



Frozen-in condition [eql] 



Po Pp 

and 



(Al) 



Polytropic relation [eq2] 

PoPo " = PpP^ ■ (A2) 

where 7 is the speciflc heat ratio. We assume 7 = 5/3 (monoatomic ideal gas). 

There are several discontinuities in the reconnection outflow, i.e., X-shaped slow shock, reverse fast 
shock, contact discontinuity, and forward V-shaped slow shock (see flgure 3). We set jump conditions 
for both sides of each discontinuity: 

X-shaped slow shock Rankine-Hugoniot (R-H) jump conditions 
Pressure-jump [eq3] 
Pi 



1 + sin 9vxp + cos 9vypY{X — 1) x { ^ — (cos OB^p + sin 9Byp)'^/2 
X [-2Vlj,^X + (- sin Ov^p + cos eVypf{X + !)]/[((- sin Ov^p + cos Ovypf - XV^p^fiipp]} (A3) 

where 



^sp 



\llpplpp > 



Vapx = V (- sin dB^p + cos OBypf j {jj,pp) , 
II is the magnetic permeability of vacuum and X is the compression ratio. 



Density-jump [eq4] 

Velocity (parallel) jump [eq5] 



^ = X (A4) 
Pp 



cos Ovi - Vq ^ (- sin Ov^p + cos Ovyp)"^ - V^^^ 
sin 9vyp - Vq (- sin Ov^p + cos Ov'^p - XVlpx) 

where Vq — {vypB^p — BypVxp)/{Bxp sin 9 — BypCOsO) is the shift speed of the de Hoffmann- Teller 

coordinate. 

Velocity (perpendicular) jump [eq6] 

-sin^., 1 ^^g^ 



— sin 9vxp + cos 9yyp X 

Magnetic held (parallel) jump [eq7] 

sin 9Bi _ [(- sin 9vxp + cos 9vypy - V2p^]X 

cos 9Bxp + sin 9Byp (— sin 9vxp + cos 9vyp)^ — XVjp^ 



(A7) 
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Magnetic field (perpendicular) jump [eq8] 

cos 9Bi 



sin 9Bxr, + cos 9By 



= 1 (A8) 



The compression ratio X is defined by the following equation (3rd order algebraic equation for X), 

{{— sinOBxp + cos 9 Byp^ [{cos 9 Bxp + sm9Byp)'^{'y — l)pp{— sm9vxp + cos9vyp)'^ 
+ (— sin 9Byp + cos 9 Bypy{2'jPp — Pp{— sm9vxp + cos9vyp)'^ + 'ypp{— sin 9vxp + cos 9vypy)]}X^ 
+{—pp{— sm9vxp + cos 9vyp)^[{— sin 9Bxp + cos9BypY['^ — 1) 
+(cos 9Bxp + sin 9Byp)'^{-f — 2)ppp{— sin 9vxp + cos 9vypY 
+ (- sin 9Bxp + cos 9BypY{{cos 9Bxp + sin 9BypY{'^ + 1) 
+2//(27Pp — pp{— sin ^i^j-j, + cos 9)vypY + ^Pp{— sin ^i^j-j, + cos 
+{//Pp(- sin 9vxp + cos 6'v^p[(cos 9Bxp + sin 9BypY^ + 2(- sin 6*53.^ + cos 9BypY{'y + 1) + 
(27Pp - pp(- sin 9vxp + cos ^%p)^ + 7Pp(- sin ^■Uajp + cos 9vypY)\}X 
+{-(7 + l)AtVp(- sin 9vxp + cos ^^2,^)^} = . 

This equation has three roots. We must choose a real positive root larger than unity. 
Reverse-fast shock R-H conditions 

Pressure-jump [eq9] 

^ = C/ (A9) 

where 

C/ = 7M1V1 - i/C/) + (1 - eJ)M + 1 

with 



and 

Density-jump [eqlO] 



= {-I + + 2M(2 - 7)(7 + l)7Mf;)/(2/A(2 - 7)) , 

Z = 7(lM + l) + (7-l)7MiV2, 
A = P^/{Bl/{2ix)) , 

Mif = (t;i - XfVA)/cf , 

Cf = \flKfpi 
Va = Bq/ y/pp^ . 



(AlO) 



Velocity jump [eqll] 
Magnetic field jump [eql2] 



Pi 

Vi — Xf 
V2-Xf 



(All) 



# = (A12) 



Si 

We must impose local magnetic flux conservation on both sides of region p and region 1. 
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Magnetic flux conservation at X-point [eql3] 

—vxpByp + VypBxp + viBi = (A13) 

Force balance at contact discontinuity [eql4] 

V-shaped forward slow shock R-H conditions 

Pressure-jump [eql5] 

P3 = Po|l + ^(a;.sin0)2(X;,-l) X 



c 



-sO 



1 {Bp cos 0)^ -2V2,,Xj, + {x, sin + 1) 

Xh 2 [(x, sin (py - XhVjoJiipo 



where c^o = \JlPo/ Po, Vaox = \l (~ sin OBq)"^ / (/xpo) and X/j is the compression ratio. 
Density-jump [eql6] 

Po 

Velocity (parallel) jump [eql7] 



(A15) 



-^Xn (A16) 



(A17) 



(A18) 



(A19) 



cos (t){v^i - XsVao) + sin ^-u^g ^ ^^ox " ^lox 
-VaqXs cos t;^ox - XuVIq^ 

where v^oa; = sin and Vao = B^j ^fjxp^. 
Velocity (perpendicular) jump [eql8] 

- sin ^{v^2, - x^Vao) + cos ^Vy-j _ J_ 
V^iXs sin 

Magnetic field (parallel) jump [eql9] 

cos (/)E^3 + sin0^y3 ^ (^^^ - Vaox)^^ 

So cos -"mOx - ^/^V^ioo: 

Magnetic field (perpendicular) jump [eq20] 

-sin0E^3 + cos (\)By'i ^ ^ 

-So sin ^ ' 

The compression ratio X/j is a solution of the following third order algebraic equation: 

{(So sin (/))^[(So sin 0)^(7 - \)pq{xs sin 0)^ + (So sin 0)^(27Po - Pq{xs sin 0)^ + 7^0(2^5 sin 0)^)]}X^ 
-F{-po(a;s sin 0)2[(So sin 0)^(7 -Fl) (So cos 0)^(7 - 2)p,pp{xs^va.(^Y 
+ (Sosin0)2((SoCos0)2(7 + 1) + (2/i(27Po - Po(a;s sin0)2 + 7po(x, sin^)^)))]}^! 
-F{//p^(a;s sin 0)^((So cos 0)^7 + 2(So sin 0)^(7 + 1) + ;Li(27Po - p^ixs sin^)^ + -fpQ{xs sin^)^))}^^ 

+{-(7 + 1)p'^pI{xs sin 0)6} = . 

We must choose a real positive root larger than unity. 
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The tip of the reconnection outflow touches the FRWF (see figure 3), hence A[ = at this point. 
This condition reduces to the following equation. 
Boundary Condition at the tip of the outflow [eq21] 

— 5^3(1 — Xc) + BxsKl — Xs) tan(f> — Xct&n6] — Bo{l — Xg) tan(j) = (A21) 

From magnetic flux conservation, the injected magnetic flux must be redistributed in the reconnection 
jet. This leads to the following equation. 
Magnetic flux conservation [eq22] 

-BypVxp + VypBj:p + [xfBi + {xc - Xf)B2] = (A22) 

We assume the following trivial relations: 

Definition from the contact discontinuity 

We can solve [eq2] for Pp, [eq8] for Bi, [eq9] for P2, [eqlO] for p2, [eql2] for B2, [eql3] for Vyp, [eql5] 
for P3, [eql6] for pa, [eql7] and [eql8] for Vx3 and Vys, [eql9] and [eq20] for B^s and By^ by hand, then 
substitute them into other nine equations ([eql], [eq3], [cq4], [cq5], [cq6], [cq7], [cqll], [eql4], [cq21] 
and [eq22]) for the following unknowns: Pp, v^p, Pi, X/, and Xs- The only parameter 

included in this problem is the plasma-/? value at the asymptotic region. By using a Newton-Raphson 
routine, with an initial guess for these unknowns, we obtain converged solutions. The procedure to 
obtain the series of converged solutions is discussed in section 2.3 in detail. 

As I Vyp\ decreases, the strength of reverse fast shock reduces, and then the pressure jump Q and 
density jump simultaneously become unity at the critical value of \vyp\. At that point, the reverse 
fast shock vanishes. We convert the coupled equations to@another set compatible to the situation with 
no fast shock, i.e., several equations reform to the following equations: 
[eq9] is replaed to [eq9'] 

[eqlO] is replaed to [eqlO'] 

P2 = Pi , 

[eqll] is replaed to [eqll'] 

^^1 = ^^2 , 

[eql2] is replaed to [eql2'] 

Bi = B2 , 

and [eq22] becomes to be equivalent to [eql3] and removed from the set of coupled equations. The 
number of coupled equations reduces to 21. This is consistent with that the locus Xf of the fast shock 
is no longer included in the set of unknowns and total number of unknowns becomes 21. 
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Fig. 1. — Spontaneous inhalation speed toward the X-slow shock. The y— component Vyp of the velocity at the region 
p versus the ambient plasma-/? for the case Vxp = (non- converging inflow) is shown. Similar to this result of Vyp, we 
can find any unknown variable as a function of /3o- These solutions for arbitrary /?o can be treated as initial guess for 
the procedure in subsection 2.4. This figure is almost the same to figure 6 of Nitta 2004, but the following points are 
different: much higher converging precision, precise treatment for junction between fast shock regime and no-fast shock 
regime. 

Fig. 2. — /3o-dependence of the reconnection rate R* . The reconnection rate is almost constant (~ 0.05) as a bmction 
of Po (plasma- /3 at the asymptotic region far outside the current sheet). The obtained value of the reconnection rate is 
consistent with previous theoretical and numerical works, but we must note that this value is spontaneously determined 
by the reconnection system itself in a frame work of ideal (non-resistive) MHD. This figure is almost the same to figure 
2 of Nitta 2004, but the following points are different: much higher converging precision, precise treatment for junction 
between fast shock regime and no-fast shock regime. 

Fig. 3. — Schematic figure of the reconnection outflow. Several discontinuities form in the outflow. We consider a 
quasi-one dimensional shock-tube problem along the x— and y— axes. 

Fig. 4. — Velocity and magnetic field at the region p versus magnetic Reynolds number This is the case for 

f3o = 0.01. These quantities are essential to determine the reconnection rate R*. As i?*^ increases, first, the transverse 
component B^p of the magnetic field quickly increases. Then, Vxp drastically developed. The negative value of Vxp 
denotes a converging inflow toward the X-slow shock. 

Fig. 5. — Direction of velocity and magnetic field at the region p. This is the case for (3o = 0.01. The circular dot and 
square denote Vyp/vxp (tangent of the inclination angle of the infiow velocity) and Byp/Bxp (tangent of the inclination 
angle of the inflow magnetic field), respectively. This graph clearly shows that the infiow velocity Vj, and the inflow 
magnetic field Bp tend to be parallel to each other when the magnetic Reynolds number increases. This causes the 
decrement of the reconnection rate ii* = Vp x Bp. 

Fig. 6. — Magnetic Reynolds number dependence of the reconnection rate R* . This is the case for /?o = 0.01. As 
increases from the case Vxp = (spontaneous inhalation), first R* increases because Bxp increases (see figure 4), then 
R* decreases as oc 1/R*^. 

Fig. 7. — Schematic picture of channel flow structure. DR: diffusion region, SS: Petschek-like slow shock. SSL: 
separatrix stream line which divides the stream lines reaching the Petschek-like slow shock and the stream lines reaching 
DR. Solid lines show the magnetic fleld lines. Arrows show the flow velocity. This is the case of significantly converging 
inflow. The right side of SSL is filled with uniform high speed converging inflow, while the left side of SSL is fllled with 
very slow speed inflow. 



